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Abstract 



We investigate the occurrence of a ferromagnetic phase transition in 
high density hadronic matter (e.g., in the interior of a neutron star). This 
^ \ could be induced by a four fermion interaction analogous to the one which 

is responsible for chiral symmetry breaking in the Nambu-Jona-Lasinio 
model, to which it is related through a Fierz transformation. Flavor SU{2) 
and flavor SU{3) quark matter are considered. A second order phase tran- 
sition is predicted at densities about 5 times the normal nuclear matter 
density, a magnetization of the order of 10^^ gauss being expected. It is 
also found that in flavor SU{3) quark matter, a first order transition from 
the so-called 2 flavor super-conducting (2SC) phase to the ferromagnetic 
phase arises. The color-flavor-locked (CFL) phase may be completely hid- 
den by the ferromagnetic phase. 

1 Introduction 

The properties of a highly compressed fermion fluid, especially with respect 
to magnetization, are of great interest since they are relevant for stellar objects 
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such as neutron stars and, possibly, quark stars [1]. Recently, magnetars, which 
are a kind of neutron stars exhibiting extremely powerful magnetic fields, have 
been discovered [2, 3]. 

We investigate a model of the magnetization of a quark fluid which is charac- 
terized by a four fermion interaction [4] related to the interaction of the standard 
Nambu-Jona-Lasino (NJL) model through a Fierz transformation. In the present 
model, a ferromagnetic transition is induced by an analogous mechanism to the 
one which is responsible for chiral symmetry breaking in the NJL model [5, 6, 7, 8]. 
This model accounts for the ferromagnetic transition at very high densities, after 
chiral symmetry restoration. This is at variance with ref.[9], where it is shown, 
in the framework of the chiral model of pions and quarks, also known as the 
linear sigma model, that spin polarization of hadronic matter may arise as a 
consequence of pion condensation, if chiral symmetry has not yet been restored. 

It should be pointed out that, in terms of a Skyrme force, a ferromagnetic 
transition has been obtained for nuclear matter [10, 11, 12]. However, in that 
case, microscopic calculations do not predict a ferromagnetic transition, at least 
below p = 7po [13]. A ferromagnetic phase of quark matter, has recently been 
investigated by several authors [14, 15, 16], who focus on its QCD origin. 

2 Symmetrical quark matter 

The model we are considering applies to high densities where chiral symmetry 
has been restored. Hence, the quark masses are reduced to their bare masses, 
which are assumed to vanish. Our model is of the NJL type [5, 6, 7, 8], but, 
for simplicity, we neglect the typical scalar-scalar and pseudoscalar-pseudoscalar 
interaction which is responsible for chiral symmetry breaking. This is because, 
at very high densities, it has no effect. Instead, we construct a pure fermionic 
QCD inspired model based on a chiral symmetric tensor-tensor interaction, so 
the Lagrangian reads 

We stress that this interaction may be understood as arising from a the Fierz 
transformation of a standard NJL interaction. We consider /i, z/ of the form 
hj ^ {1; 2, 3}, i 7^ j. In the following we will consider massless quarks and we will 
work in the no-sea and mean-field approximations. By mean-field approximation 
we mean that the operator V'7i7iTfeV' is replaced by {tp'ji'jjTkip) +Tp'~ii'~ijTk'ip — 
{'4'li'lj'^ki^) and second order terms (and higher) in ip'ji'yjTkip — {^lilj'Tki^) are 
neglect. In the Conclusions, the implications of the no-sea approximation will be 
briefly discussed. Let us assume polarization along the 3 axis. Then, (^7i72Tfc'?/')^ 
is replaced by 2{%l}'-^i'-^2Tki^)''Plil2'Tki^ — {'4'lil2'rk4')'^ ■ We introduce the notation 

Fk = Gi^E^Tk^) = 2G'(^7i72rfc^). (1) 



Manifestation of this interaction is seen only at very high densities. In the mean 
field approximation, the postulated interaction leads to the Lagrangian density 

C-MFA = ^(Jpr^,^P) - Fki^E.TktP) - ^. (2) 

In turn, this leads to the Dirac equation 

{-ia ■ V + e^F/3E3)V' = eijj, 

where e,- = 1 for quarks u and e,- = — 1 for quarks d denote the eigenvalues of Ts, 
and Fk is a parameter which is determined from (1). In the following, we take 
Fi = F2 = 0, and F3 = F. 

In this section we will consider symmetric 2-flavor quark matter, the Fermi 
energy being the same for quarks u and d. The single particle energy eigenvalues 
satisfy 




e, = ±J{\F\±Jpl + pl] +pI, (3) 



where the ± signs are such that the single particle energies of quarks u with 
momentum p and spin s are the same as the energies of quarks d with the same 
momentum and spin —s. By ^3 we denote the component of momentum along 
the direction of spin polarization, pi and p2 being the transverse components. A 
no sea approximation shall be used. Thus, the quark Fermi surfaces are given by 



f^=[\F\±^Pi + pij +Pi, 

the + sign corresponding to the Fermi surface of quarks u with spin up or quarks 
d with spin down, (in the direction of magnetization) and the — sign to the 
Fermi surface of quarks u with spin down or quarks d with spin up. Since the 
temperature is assumed to vanish, the Fermi energy coincides with the chemical 
potential fi. 

We observe that the formalism behind the description of ferromagnetic con- 
densation is analogous to the BCS theory (the parameter F playing the role of the 
gap A), or to the NJL theory for chiral symmetry breaking (the order parameter 
F/G corresponding then to the quark condensate (qq))- The condition 

replaces the gap equation. Below the critical potential Hc, this condition implies 
that F = 0, so that the state of equilibrium is the normal phase. Above Hc, 
equilibrium occurs for F ^ (cf. Fig 1). The critical chemical potential /Xc is 
determined by the condition 

d'Hfi,F) ^ 

dF^ F=o 
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Figure 1: The thermodynamical potential $ vs. F, for symmetric flavor SU{2) 
quark matter, for fixed values of //, being /ic = 0.41 GeV. 



and is given by 



A^c 



TT 



3G 



The Fermi surface may be visualized as a circle of radius fi rotating around a 
line lying on its plane at the distance F from its center. If the line lies outside the 
circle, the Fermi surface is doughnut-shaped and we have full polarization. If the 
line crosses the circle but is not a diameter, we have two Fermi surfaces, one for 
spin up and the other one for spin down, corresponding to partial polarization. 
If the line is a diameter, we have no polarization at all, since the Fermi surface 
is the same for particles with spin up or down. Thus, two different expressions 
are obtained for the thermodynamical potential $(F, /i), one which is valid for 
\F\ > fi (needed at higher densities), and another one for |F| < fi (needed at 
lower densities). 

Taking into account that the color- flavor degeneracy is 6, and keeping in mind 
that 27r^/i^F is the volume of the torus generated by a circle of radius /i whose 
center describes a circle of radius F, it is easily found that, for \F\ > fi, when 
full polarization is present, the particle number reads. 



N = 

and the energy reads 
QV 



6V 



|dV(^,-^) = ^vrVl/^^ 



(5) 



E 



(2vr) 



d^p9{ep - fi)ep + 



VF^ 



2G (27r 



8V 2,r., s VF^ 



(6) 
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Figure 2: The order parameter F/G versus the chemical potential /x, for sym- 
metrical matter (SM). In the normal phase, F/G=0. 



where the term VF'^/{2G) comes from the standard mean field approximation. 
The limits of integration over the transverse momentum ^/pf+p^ in Eqs. (5) 
and (6) are F — fi and F + fi. The thermodynamical potential becomes 



^ = E-Nfi 






(27r 



This is the expression which is valid for F > n. The physical F value minimizes 
this expression and is equal to 



4G 



(2^ 



so that in the fully polarized phase the thermodynamical potential reads 



$ = - 



V G/i« 



7r2 8 



However, this expression is valid only for f^ > fj^p = \/2n/G, since \/h\jG is the 
\x value such that minimum of $ occurs precisely for F = [i. For /i < [i^^ the 
condition (4) implies that F = 0, so that the state of equilibrium is the normal 
phase for which the thermodynamical potential equals 



$ = - 
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Figure 3: The pressure P vs. the chemical potential /i for G = 20 GeV~^ 
comparing the ferromagnetic phase (FP) and the normal phase (NP) of symmetric 
matter. 

In order to describe the system for \ic< \i < /Xp, we need $(/x, F) for < F < 
/i. This is the case of partial polarization, which we now briefly discuss. Then, 
the limits of integration over the transverse momentum \/v\ + v\ ci-re 0, fx — F 
(spin up) and 0,fi + F (spin down), so that the quark number equals, 

Ni^^, F) = -^, (y"d3p^(^P - /i) + j^'pO{e, - /i) 

= ^ I V/i^ - F\F^ + 2/i2) + 3F/x2 tan-i , ^ ) . (7) 



TX 



^fJF^^F^ 



For the energy we find 

^ ^ VF 
E = El 



2G 



where 



Eiifi, F) = — -^ ( / d=^p Ep e{ep -^i)+ d=^p Ep e{ep - /i) 



d/i /i^:— = ^^N — djj, N 
9fi Jo 



= fiN- dfiN- d/iAr = /iAr--^F^ + $i, (8) 

Jo Jf °^ 

In the case < |F| < yU, the thermodynamical potential $ is given by 

41/ VF'^ 

*(,.f) = *,--f^ + — . (9) 
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Figure 4: The pressure P vs. the quark density p for G = 20 GeV~^, for symmet- 
rical quark matter, in the ferromagnetic (FP) and in the normal (NP) phases. 
Partial polarization is described by the portion of the graph with lowest slope, 
which occurs for 4.87po <P< 9.99po and for 0.00175 < P < 0.00273 (GeV^). 



where. 



/^ 



VF^ /2p3 + 3^F2 



*.(,. F) = - y^ d, ^ = -^ ^-j_^ y^iT^ 



,3 



1 p+v//i2-F2 



|F|3 ^^2 _ p2 A ^ \F\ 2) ^ ' 

In (9), \F\ is determined by minimizing $ for fixed /x. It is found that for /x < iic, 
the minimum of $(p, F) occurs for F = 0, so that, for low densities, the normal 
phase prevails. For p^ < p < ^J2t^ /G = ^p the minimum of $(p, F) occurs for 
< \F\ < II, and then partial polarization is realized. 

Numerical results for G = 20 GeV~^, arbitrarily chosen but of the order of 
the coupling constant of the chirally symmetric four-fermion interaction of the 
NJL model which takes values ~ 5 — 10 GeV~^ [8], are presented in Figs. 2 to 4. 

In Fig. 2, the order parameter F/G is plotted versus the chemical potential 
p, for symmetrical matter. The ferromagnetic phase occurs for a quark chemical 
potential above fic =0.41 GeV. A transition of second order, characterized by a 
continuous behavior of the magnetization and of the baryonic density, is found. 

In Fig. 3, the pressure P is plotted versus the chemical potential p for both 
the ferromagnetic and the normal phases of symmetric quark matter. The most 
stable phase is the one having the largest pressure for a given chemical potential. 
Therefore, for the lowest chemical potentials the normal phase is the equilibrium 
phase for both symmetric and neutral matter, while at high chemical potential 
the ferromagnetic phase is the stable one. 
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Figure 5: The critical values of the pressure and of the chemical potentials versus 
the couphng constant G. By ftp we denote the critical chemical potential for the 
onset of full polarization. 

In Fig. 4, the pressure P is plotted versus p/po, where p is the baryonic density 
and po is the nuclear saturation density, comparing the ferromagnetic and the 
normal phases of symmetrical quark matter. Partial polarization is described 
by the portion of the graph with lowest slope. For symmetrical matter, this 
region occurs for 4.87po <P< 9.99po, and 0.00175 < P < 0.00273 (GeV^). The 
critical quark densities are quite high. A larger coupling constant would induce a 
transition at lower densities. In Fig. 9, the dependency of the critical values of the 
pressure and of the chemical potential on the coupling constant G is displayed. 
A phase transition at densities comparable to the ones expected to exist in the 
interior of neutron stars, twice to four times the equilibrium density of nuclear 
matter, would require that a G of the order 25 GeV~^ to 30 GeV~^. Although 
these couplings are probably too large, they may be justified, as explained below, 
if the vacuum polarization is taken into account. Only values of p up to lOpo are 
relevant. 



3 Magnetization 

As suggested in [4], magnetization of a relativistic fermion fluid may be asso- 
ciated with the lowest Landau level (LLL). The energy of the Landau level I'jPs is 
obtained through the replacement pl+ P2 ^ 2v\Q\B in the expression (3) of the 
single particle energy, so that, in the presence of a magnetic field 5, the energy 
of a quark with charge Q is equal to 



eu{pz 



[f±^2\Q\Bu 



Vi- 



la our model, the effective LLL does not arise for i/ = 0, but for v = F'^ / {2\Q\B) , 
so that its energy is eiiLiPs) = bsl- By analogy with [4], we postulate, some- 
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what ad hoc, that the contribution of the effective LLL to the thermodynamical 
potential reads 



p-i 



where \Qr\ is the quark charge and \Qr\B /{2ti) the Landau degeneracy. The 
magnetization of fully polarized matter is then equal to 

, . _ ^^^\\Qu\ + ig^l) _ 3/i^(||e| + ||e|) 



in^ 



87r2 



For fic < fJ' < f^p ^e will have 
M 



In fig. 6, the magnetization is plotted vs. the chemical potential, for symmetric 
quark matter. 



4 Neutral quark matter 

In electrically neutral quark matter, the density pd of quarks d is twice the 
density pu of quarks u. In the normal phase \F\ = and the numbers of quarks 



u and d are given by 

V . V . 

Nu = —K. Nd = —fii 

where /i„, fid are the Fermi energies of quarks u, d, respectively. Since pd = 2p„, 
we have pd = 23/i„. The thermodynamical potential is written $ = K^ + Kd — 
fi{Nu + Nd). Here fi is the chemical potential which fixes the quark number 
independently of fiavor and the kinetic energies are equal to Ku = ^/^t, and 



Kd = 4^/irf, so that 



* = 4^^/"- + ^'d) - /^^(/^« + /^d)- 



The values of /i„, pd are not fixed a priori but are determined by minimization of 
$ under the condition /i^ = 23 fi^, which leads to 

_ 3 _ 3-23 

f^u — r/^! f'd — -fJ^, 

1 + 23 1 + 23 

so that 

4vr2 V1 + 2I; 

In the fully polarized ferromagnetic phase, |F| > pd (since fid > fiu) we have 
(cf. eq. (5)) 

N^ = —\F\fil Nd = —\F\fil 
Therefore, for electrically neutral matter we obtain fid = \/2fiu, and 

K^= \F\fil Kd = -\F\fil 
The thermodynamical potential reads 

'^ = ^\FM + f^l)-f^^\F\{fil + f^l) + ^. 

Under minimization with respect to fiu-, fJ^d subject to the condition fid = \/2fiu, 

we find 

3 

f^u — ~ ~tA* 

1 + 22 

so that 

47r(l + 2v^)2'^ 2G ^ ^ 

Upon minimization of the previous expression with respect to F we obtain 

25(l + 2v^)%2^ ■ 
10 



In order to describe the regime of partial polarization, the condition that for 
each quark u there are two quarks d is imposed by requiring that 

^N{fi^,F) = N{fid,F), (13) 

where the function N{fi,F) is given by (7). The thermodynamical potential of 
neutral quark matter equals 

$(/i, F) = - (Ei(/i„, F) + Ei(/i,, F)) + ^^ - /i (^- {N{^.,, F) + iV(^,, F)) 

(14) 
where the functions A^(/i, -F), Fi(//, F) are given by (7), (8). The factor \ accounts 
for the flavor degeneracy which was included in the definitions of these functions. 
The equilibrium thermodynamical potential is obtained by minimizing (14), for 
fixed /i, with respect to /i„, /i^, and F, under the constraint (13). 

We find that the behaviors of neutral quark matter and symmetric quark 
matter are qualitatively very similar, the pressure of neutral quark matter being 
slightly smaller, for a given /i, due to the charge neutrality constraint. 



5 Strange quark matter 

In this Section, we consider strange quark matter at densities which are so 
high that chiral symmetry has been restored, in the sense that the masses of 
quarks m, d and s are reduced to the bare masses. Moreover, it is assumed for 
simplicity that the bare masses vanish. For quarks s, this is admittedly a very 
drastic assumption. However, we will analyze its consequences and comment at 
the end. As in Section 2, we postulate a QCD inspired tensor-tensor interaction, 
so that the Lagrangian reads 

£ = i^^^d^i, - ^G'(V^7'^7"AfcV')(V'7A.7.A,V^), 

where A^ denote the S't/(3) flavor Gell-Mann matrices. This Lagrangian is not 
chiral invariant. The terms which insure chiral invariance have been omitted, 
for simplicity. We assume polarization along the 3 axis, consider the mean-field 
approximation, and define 

Ffc = iG{^-ixl2^ki>)- 

In the mean field approximation, for k ^ 3, 8, we have Fk = 0, so that the 
Lagrangian density reduces to 

CMFA = t{i^rd,i^)- Yl Fk^sAk^p)- Yl §• 

fcG{3,8} fcG{3,8} 
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In turn, this leads to the Dirac equation 

{-idt + iajVj)ij-l3T.^ J2 FkKk^ = Q. 

A;6{3,8} 

The operator 'Yl,k<^{z 8} ^k^k has different values for quarks m, c?, and s. Its eigen- 
values are equal to 

J-. = (f^ + -^Fg") 5r,u - (fs - -^Fg") 6r,d - 4^Fs6r,s (15) 

so that the Dirac equations reduce to 

{—ia. ■ V + Jv/JSs)?/',- = Er-ipT, T G {u, d, s}. 

We assume that J^u is positive and J-^, J-'g are negative. The single particle energy- 
eigenvalues read 




^p.^ = M [\^r\± ypi + pi] +pi 



n2 

5 



where the + sign refers to quarks u with momentum p and spin up and to quarks 
d, s with the same momentum and spin down, while the — sign refers to quarks u 
with momentum p and spin down and to quarks d, s with the same momentum 
and spin up. The quark Fermi surfaces are given by 



fi'= l\Tr\±^Jpi+P2j +P3 , 

the + sign corresponding to the Fermi surface of quarks u with spin up or quarks 
d, s with spin down, and the — sign to the Fermi surface of quarks u with spin 
down or quarks d, s with spin up. 

Several different analytic expressions, each one with its own range of validity, 
are obtained for the thermodynamical potential. One which is valid when |J-V| > 
^, T G {u,d,s}, (needed at high densities), and another for |Jv| < /i (needed at 
low densities). Another one which is valid when the largest |J-V| is bigger than n 
(needed at high densities), and the remaining |J-V|'s are smaller than /i. Here, fj, 
is the Fermi energy, which is identified with the chemical potential. Taking into 
account that the color degeneracy is 3, in the first case, when full polarization is 
present, the r-quark number reads, 

^- = (Sy^/^'P ^(^P,^-/^) = ^3^"\^r\^^' , (16) 



and the r-quark kinetic energy reads, 

3V f 

^r = 77^3^ / d^P 0i£p,T - /i) £p,T = t^zT^'^'^IJ^tIi^'^ ■ (17) 
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Figure 7: $ vs, F3, for strange quark matter. The curvature at the origin turns 
upside-down for fi about 0.410 GeV. 

The total quark number and total energy read, respectively 

QV 



N 



(2vr)= 



■IT 



,2,2 



/i" (|J^„| + \Td\ + \J^s\), 



E = ^vrVd J-„| + I J-,1 + I J-,1) + ^(F| + Ff 



s;) 



where Jv is given by (15) and the term V{F^ + F^)/{2G) comes from the stan- 
dard mean field approximation. For full polarization, the thermodynamical po- 
tential becomes 



^ = E- Nfi 



^vr^d^J + l^.l + l-F.iy + ^{F^ + F^) 



As a simplifying ansatz, let us assume that 



This assumption places the quarks d and s on an equal footing. Then, 



*' n Q-'^Sy '' d •'s 



3 ^ 



(19) 
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Figure 8: P vs fi for strange matter. The crossing point of the 2SC and spin 
polarized phases is lower for the higher G. 



Under this ansatz, for ^F^, > ji (which insures full polarization, i.e., | J-'„|, \J^d\i \^s\ > 
/i) the thermodynamical potential equals 

161^ 



$ 



3(27r)3 



2V 
3G^' ' 



implying that 



n'\F,\^^' + ^F^ 



,f^^ 



(20) 



Fg = G^ , 

^ 27r ' 



so that in the ferromagnetic phase, which is realized if fi is above the critical 
value fip, the thermodynamical potential reads 



$ 



V 



67r2 



G/i^ 



(21) 



This should be compared with the expression of the thermodynamical potential 
in the normal phase which is realized if /i is below the critical value ^c and equals 



$ 



3V 



/. 



(22) 



Then, in the regime of partial polarization, the thermodynamical potential of 
strange quark matter equals 



^(/i, F3, Fg) = -($i(/i, J-„) + $i(/i, J-,) + $i(/i, J-,)) 



(23) 
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Figure 9: The pressure versus the baryonic density, comparing the 2SC phase 
with the spin-polarized phase. 

where $i(//, F) is given by (10). Under the ansatz (18), (19), for < ^F^ < /i 
(or < 1^3 < |/i), the thermodynamical potential reduces to 



$(/i, Fs) = i Ui L, ^\Fs\] + 2$i L, ^iFsl 



8V 
97r 






(24) 



For 1^3 < /i < 1^3 (or |/i < |F3 < /i), the thermodynamical potential reduces 
to 



2 \ 81/ 

$(/X, F3) = $1 (/i,-|F3|j -^1^3 



1^'^' + ^§- '^^) 



^J4 " i;. ,,3 , T/-^3 



Therefore, we are faced with three partial expressions. Of course, the equilibrium 
thermodynamical potential is obtained by optimizing these partial expressions in 
their respective domains of validity. We find that jJc = vt/a/SG is the critical fx 
for the spin polarization phase transition, and ^p^ = ^^/3^^/G is the critical /x for 
the onset of absolute full polarization. The expression of the critical chemical 
potential /ip^ for the onset of full polarization of quarks u (but not of quarks d, s) 
is not given because it is cumbersome. 

Color superconductivity in quark matter has been investigated in [17], where 
a color-flavor symmetric version of the so-called 2 flavor super-conducting (2SC) 
phase has been proposed and its transition to so-called color-flavor locked phase 
(CFL) was analyzed. In [17], color superconductivity has been described by a 
mean-fleld thermodynamical potential operator (i.e., the Hamiltonian constrained 
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Figure 10: The minimizing F^ vs /i, for strange quark matter, comparing G = 20 
GeV~^ with G = 25 GeV~^ The critical /i is lower for the higher G. 

by the quark number conservation) of the form, 

Km FA = ^^pcla 



^iap^iap 



lap 



~r / ^ \ kc jbp ^iap ^ijk^abcT]p ~r ^kc ^iap ^jbp ^ijk^abc'lp) 



^kc 



ijkabcp 

'TT / J \^jbp Ciap) ^ijk^^abcVp^ 
ijabp 



V 
2G, 



jci 



(26) 



jc 



where e^ = a/pM-M^— /i, M being the current mass. In the 2SC phase presented 
in [17], the gaps A^^ are real and independent of the color-flavor indices, that is, 
Ajc = A, j G {M,(i, s},c G {r, (yf,6}, so that the two specific flavors or the two 
specific colors which are paired are not pure colors or flavors, but are related 
to pure colors and flavors through a certain rotation in color-flavor space, which 
has the great advantage of insuring automatically color-flavor symmetry. In the 
CFL phase, the indices j, c are coupled in a prescribed way, for instance, {u with 
r), [d with g), [s with h) (cf. [17] for further details). For the model considered 
in [17], the 2SC phase is stable for yU < 0.505 GeV, while the CFL phase is stable 
for /i > 0.505 GeV, for Gc = 5GeV~^. In that model, the quarks have all the 
same helicity (either positive or negative). If we wish to take into account both 
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helicities, the gaps reported in [17] are valid and the pressure becomes twice the 
pressure reported there, under the condition that half the value of Gc is used, 
that is, for Gc = 2.5GeV-^. 

It is of interest to investigate the interplay between color superconductivity 
and spin polarization as is here described. The transition between the 2SC phase 
and the ferromagnetic phase is depicted in Fig. 8, where the pressure P is plotted 
vs the chemical potential fi for the normal phase, for the ferromagnetic phase, 
for the 2SC phase and for the CFL phase. The curves corresponding to the spin 
polarized F phase, to the 2CS phase and to the CFL phase seem to cross at the 
same point. This is by pure accident. If instead of G = 20 GeV~^ a slightly 
larger value is considered, the F curve will cross the 2CS and the CFL curves 
before these curves cross each other. This is clearly seen in fig 8 (b), where the 
spin polarized phase for G = 25 GeV~^ is plotted. In our model, the CFL phase 
seems to be completely hidden by the ferromagnetic phase, but, if a slightly lower 
G is used, the CFL phase will be realized for a short while. Spin polarization 
and color superconductivity has recently been considered in [18]. The relevance 
of the interplay between the CFL phase and the 2SC phase for the stability of 
hybrid stars is discussed in [19]. 

Suppose we have an unperturbed Lagrangian Cq and two types of perturba- 
tions, Cj and Cji, so that Cq + Cj describes the phase / associated with some 
condensate and £o + '^ii describes the phase II associated with another con- 
densate. The interplay between the two phases is described by the Lagrangian 
Cq + Cj + Cij. However, in order to treat approximately this interplay we may ap- 
ply the Gibbs criterium to the first order transition between the phases described 
by the truncated Lagrangians. Clearly, it would be more correct to use the full 
Lagrangian, following [20], and then the first order transition may be replaced 
by a second order one. We have used the approximate method to describe the 
interplay between the 2SC phase and the spin polarized phase. 

6 Conclusions 

Spin polarization of a high density hadronic fluid of quarks has been studied. 
We describe a spin polarized phase in hadronic matter arising from a four Fermion 
interaction which may be regarded as the relativistic analogue of the interaction 
considered in [10, 11, 12] and is related to the interaction of the standard NJL 
model through a Fierz transformation. A no sea approximation has been used. 
Flavor SU{2) and flavor SU{?>) quark matter are investigated. A second order 
phase transition at densities about 5 times the normal nuclear matter density is 
predicted and a magnetization of the order of 10^^ gauss is expected. Numerical 
results for G = 20 GeV~^, which is arbitrary but of the order of the coupling 
constant in the NJL model, predict a phase transition at 5—8 times the saturation 
density po, which is probably too high. Increasing the coupling to G = 30 GeV~^ 
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would decrease the transition density to 3po- It is also found that in flavor 
SU{3) quark matter, a first order transition from the 2SC phase [17] to the spin 
polarized phase arises. The spin polarized phase is stable for fi > 0.51 GeV. 
The 2SC phase is stable for /x < 0.51 GeV. In our model, the CFL phase [17] is 
completely overshadowed by the spin polarized phase. 

We have taken G = 20 GeV~^. The critical chemical potential for symmetric 
quark matter and for strange matter turns out to be fic =0.41 GeV. If lower G 
values are considered, fic increases and in that case the transition from the 2SC 
to the CFL phase occurs before the transition to the spin polarized phase, so 
that the CFL phase is realized. If higher G values are considered, fie decreases, 
so that the transition to the spin polarized phase occurs before the transition 
from the 2SC phase to the CFL phase, which is then completely overshadowed. 
Due to the vacuum polarization, it is still meaningful to consider slightly higher 
G values. Indeed, the main effect of taking the vacuum polarization into account, 
is the appearance, in the expression of the thermodynamical potential, of a new 
term of the form — kA^F^, where k > is a certain numerical factor and A is 
the regularization cutoff. But this will only renormalize the coupling constant 
G, so that we may disregard the vacuum polarization provided we replace the 
unrenormalized Go by a new G such that — kA^ + l/(2Go) = 1/(2G), implying 
that Go < G. This may justify the use of a large G. For simplicity, we have 
assumed that the strange quark current mass vanishes. This is an acceptable 
assumption for /i^ =0.41 GeV. 

This has been an exploratory work where the possibility of a phase transition 
to a magnetized phase in quark matter was investigated. The application to a 
physical system like a compact star may require the careful consideration of the 
vacuum contribution and of finite quark current masses. 
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